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Structure of eigenstates in a periodic quasi- ID waveguide with a rough surface is studied both 
analytically and numerically. We have found a large number of "regular" eigenstates for any high 
energy. They result in a very slow convergence to the classical limit in which the eigenstates are 
expected to be completely ergodic. As a consequence, localization properties of eigenstates originated 
from unperturbed transverse channels with low indexes, are strongly localized (delocalized) in the 
momentum (coordinate) representation. These eigenstates were found to have a quite unexpeted 
form that manifests a kind of "repulsion" from the rough surface. Our results indicate that standard 
statistical approaches for ballistic localization in such waveguides seem to be unappropriatc. 
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In the past decade much attention has been paid to sta- 
tistical properties of eigenstates of closed disordered sys- 
tems. As a result, to date practically everything is known 
for quasi-lD systems with "bulk disorder". The success 
is mainly related to the developments of the non-linear 
a model (see, e.g. [1] and references therein). One of the 
most important results is that the statistical properties 
of such systems are essentially determined by one char- 
acteristic length only, known as the localization length of 
eigenstates. This fact is entirely due to a strong mixing 
between transverse channels resulting from bulk scatter- 
ing and leading to a diffusive character of transport. 

A much more difficult situation was found to occur for 
the models with surface scattering when the disorder is 
due to a surface roughness. In quasi-lD geometry such 
models are closely related to optical/microwave waveg- 
uides and have many physical applications in different 
fields [2]. The main problem in the rigorous treatment 
of this kind of systems is in the ballistic character of 
the scattering which has weaker statistical properties in 
comparison with diffusive scattering. Progress in this di- 
rection is related to recent developments of the " ballistic 
sigma-model" , however, the problem is still far from be- 
ing solved (see discussion and references in [3]). 

As was recently shown in a number of numerical stud- 
ies [4,5], the transport in quasi-lD waveguides with rough 
surfaces is highly non-isotropic in channel space. Specif- 
ically, the transport through such waveguides strongly 
depends on the incident angle of incoming wave. In par- 
ticular, the transmission coefficient smoothly decreases 
with an increase of the angle, since characterictic lengths 
for backscattering are different for different channels [5] . 

To understand generic features of surface disordered 
quasi-lD systems, in this Letter we perform a detailed 
study of the structure of eigenstates of a 2D quantum bil- 
liard (or waveguide). We consider billiards which are pe- 
riodic in the longitudinal coordinate x, and with Dirich- 
let boundary conditions on the low, y = 0, and up- 



per, y = d + a£(x), surfaces with £(z + b) = and 
(£(x)) = 0. Here the angular brackets stand for the aver- 
age over one period b (or, in the case of a random profile, 
over different realizations of £(x)). 

Our main interest is in the study of the structure of 
eigenstates of this billiard in dependence on their energy 
and model parameters. For this we use the technique that 
transforms the Hamiltonian H = 2^-{Px + Py) f° r a f ree 
particle inside the billiard with the above boundaries, to 
a new Hamiltonian which incorporates surface scattering 
effects into effective interaction potential. This can be 
achieved by the transformation to new canonical coordi- 



nates, u — x; v 



with e = a/d. As a result, the 



l+e£(s) 

boundary conditions for new wave function are trivial: 
$(u, v) — at v — and v = d (see details in [6]) 
In the new variables the Hamiltonian gets [6] , 
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and £ u = d^/du, £ uu = d 2 £,/du 2 . 

One can write the Hamiltonian in the following form, 

H = H° + V(u,v,P u ,P v ); H° = -?-(P 2 + P 2 ) (3) 

Zm e 

where P u and P v are the new canonical momenta. In 
this way, the "unperturbed" Hamiltonian H° descibes 
free motion of two "particles" inside a billiard with flat 
boundaries, y = 0, y = d, and V stands for an effective 
"interaction" between the "particles" . Such a represen- 
tation turns out to be very convenient for the study of 
chaotic properties of our model, since one can use the 
tools and concepts developed in the theory of interacting 
particles (see [7] and references therein). 
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This model has been thoroughly studied in Rcfs. 
[6,8-12] for the specific case £(x) = cos(2nx/b). The 
main interest was in the properties of energy spectrum 
[6,10], and in the quantum-classical correspondence for 
the shape of eigenfunctions (SEF) and local density of 
states (LDOS) [8,11]. In particular, it was shown that 
for highly excited states the global properties of the SEF 
and LDOS in the quantum model (3) are similar to those 
described by the equations of motion for a classical par- 
ticle moving inside the billiard. On the other hand, quite 
strong quantum effects have been revealed for individual 
eigenstates in a deep semiclassical region [11]. 

Below we address the case of a rough surface £(x) = 
J2k=i cos(2nkx/b), focusing on the properties of 
eigenstates. The surface is modeled by a large sum of har- 
monics with randomly distributed amplitudes Ah. With 
an increase of Nt the degree of complexity of £(x) in- 
creases and for a large Nt ^ 100 the surface can be 
treated as the random one. 

Since the Hamiltonian (3) is periodic in u = x, 
the eigenstates arc Bloch states and the solution of 
the Schrodingcr equation can be written in the form 
iPe(u,v) = cxp(ix) ip x {u,v), with ip x (u + 2ir/b, v) — 
tp x (u, v). Here the Bloch wave vector x(-E') is in the first 
Brillouin band, (— f < X < f )• By expanding ip x (u,v) 
in the basis of H , the a-th eigenstate of energy E a (x) 
can be written as 



oo oo 



r(«,«;x)=E £ C a mn {xW mn {u,v), 



(4) 



7r 1 l 2 g 1 ^ 4 sin(m7rv/d) exp{i(x 



where ^ n (u,v) 
2nn/b)u}. The factor g = [1 + e£( u )] 2 arises from the 
orthonormality condition in the curvilinear coordinates 
(u, v) (see details in [6]). 

In the "unperturbed" basis defined by H°, the matrix 
elements of the "interaction" V(u, v, P u , P v ) can be writ- 
ten explicitly for any profile £(x) [11]. This fact is very 
useful in the study of the dependence of the properties 
of eigenstates on the form of profile. 

The eigenvalues of H° are given by the expression, 
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In numerical simulations we have to make a cutoff for the 
values of m and n in the expansion (4) . Our main results 
refer to the ranges, 1 < m < M max and |n| < N max with 
N m ax = 32 , M max = 62, for which the total size of the 
Hamiltonian matrix is L = (2N max + l)M max = 4030. 
Since the statistical properties of eigenstates do not de- 
pend on a specific value of the Bloch index x inside the 
band [10], we fix it to x = 0.l7r/6. 

One natural representation of the Hamiltonian matrix 
Hl,l'{x) — < I I H \ V > x is the so-called "channel rep- 
resentation" for which one fixes the values of n starting 



from the lowest one, n — —N maxi with the run over m for 
each n. In this way the matrix has a block structure that 
manifests peculiarities of the interaction between differ- 
ent channels specified by the index m. For our purpose, 
however, to analyze the properties of the eigenfunctions 
it is more convenient to use the "energy representation" 
according to which the unperturbed basis is ordered in 
increasing energy, Ef +1 > [11]. 

In what follows we mainly discuss periodic billiards 
with a weak roughness, e = a/d -C 1. However, all ma- 
trix elements of the Hamiltonian are computed accord- 
ing to exact analytical expressions. This is important 
because the contribution of the "gradient" terms (that 
depend on £„, £ uu ) is strong for N -c 1 and should be 
treated non-perturbatively. 

In order to characterize quantitatively the structure of 
eigenfunctions we compute the entropy localization length 
Ih, given by the expression, 

l H = exp {- (H - Hgos)} « 2.08 exp{-W}. (6) 

Here Ti = y]f—i wf Inwf stands for the Shannon en- 
tropy of an eigenstate in a given basis, and TLgos is 
the entropy of a completely chaotic state characterized 
by gaussian fluctuations of the components C" with the 
variance < wf >= \Cf\ 2 = l/L [13]. Note, that the 
value of Ih is proportional to the localization length li pr , 
defined via the inverse participation ratio, U pr = 'SP^ 1 
with V = J2f=i( w ? ) 2 [13]- Both quantities give an es- 
timate of the effective number of components in exact 
(perturbed) eigenstates. 
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FIG. 1. Localization measure Ih for a rough surface with 
Nt — 100 for e = 0.06 and d = b; (a) for exact eigenstates 
| a >, (b) for individual LDOS | I >. 

In Fig. 1(a) the value of Ih is plotted versus the in- 
dex a for exact eigenstates ordered in energy E a . This 
typical dependence of Ih on a is quite instructive. As 
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one can anticipate, the number of principal components 
in the eigenstates increases, in average, with energy. On 
the other hand, for any large energy E a there are many 
eigenstates that have small values of Ih- For understand- 
ing the origin of these strongly localized eigenstates (in 
the unperturbed basis | I >), it is convenient to consider 
the so-called individual LDOS [8]. This quantity corre- 
sponds to the representation of an unperturbed state | I > 
in the basis of exact states \a >. Using the definition (6) 



with 7i 



1 lnui" where the sum now runs over 
a for a specific value of Z, one can characterize how many 
exact states contains specific unperturbed state | I >. 

The data of Fig. 1 (b) show that there is a large number 
of unperturbed states that seem to be close to the exact 
ones (with Ih ~ 1). The important point is that these 
states appear in a regular way as a function of I. The 
analysis shows that exact eigenstates | a > with smallest 
values of Ih are originated from the unperturbed states 
with m = 1, see Eq.(5). The data of Fig. 1(a) man- 
ifest that localized eigenstates | a > with small values 
of Ih can be classified by groups that are characterized 
by the values of m = 1, 2, 3 . . . (for small m) of those 
unperturbed states | I > to which they are " close" . 
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FIG. 2. Examples of exact eigenstates that are strongly 
localized in energy representation; (a) a — 968 (m = 1), 
(b) a = 976 (m = 2). The probability \t/j a (x, y) \ 2 is plot- 
ted for the parameters of Fig. 1 (dark regions correspond to 
high probablity). Broken horizontal curves show upper sur- 
face profiles. 



the explanation of this phenomenon by making use of the 
Hamiltonian in variables (u, v). 

We start with the fact that eigenstates of the type 
shown in Fig. 2(a) are exponentially localized in the 
m — space (channel space) independently on n. There- 
fore, one can write, C"(fc) = C^„(fc) <~ exp{— fi(m — 1)} 
where f3 is some constant (of the order of unity) deter- 
mined by numerical data, and surprisingly independent 
of energy. Therefore, these localized eigenstates in the 
(u, v) variables have the following form, 



,,■ - C p i(x+2™ a /b)u \ " 
rloc — 1/4 e l^i 
• y m=l 



e -/J(m-i) sin 



V d 



(7) 



Here C is the normalization constant determined by 
the orthonormality condition in curvilinear coordinates, 

IT !o dudv ^To*^ioc = 1 ( scc I 6 ! for details). As a 
result, one obtains, 

I , h a M |2 = e 2/3 - 1 sin 2 (7n//rf) 
Vioc\y) 4?r cosh/ 3-cos(7ry/d) y > 
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FIG. 3. Projection of the eigenfuction profile of Fig. 2(a) 
onto the y— coordinate (dots), together with the analytical 
expression (8) (white curve). The value (3 — 0.53 was numer- 
ically found by fitting C™„(fc) to the exponential dependence 
exp{— P(m — 1)}. In the inset, profiles (8) for strongly local- 
ized eigenstates with m = 1 in dependence on the number Nt 
of harmonics are shown for (a) Nt = 1, f) — » oo, (b) Nt = 50, 
P = 0.77, and (c) iV T -» oo, -f 0.48. 



Two examples of such eigenstates are given in Fig. 2 
in the coordinate representation (x,y). It is quite un- 
expected that these eigenstates are very different from 
the unperturbed ones even though they are quite close 
in energy representation. Fig. 2(a) shows that the rough 
boundary "pushes" the probability \ip a (x, y)\ 2 away from 
it, differing from the unperturbed mode with m = 1 
whose maximum is at the center, y = d/2. Similar re- 
pulsion occurs for m = 2, see Fig. 2(b). Below we give 



The comparison of this expression with the numerical 
data is shown in Fig. 3. One can see that in spite of 
a relatively weak coupling of low channels (m ~ 1) to 
all others, the scattering from a rough surface strongly 
modifies the unperturbed states in y— direction. One can 
speak about a kind of "repulsion" of such eigenstates 
from the rough surface. In order to see how this repul- 
sion depends on the degree of roughness, we have studied 
the form of the states with m = 1 in dependence on the 
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number Nt of harmonics in the surface profile. The re- 
sults shown in the inset, reveal that with an increase of 
Nt the position of the maximum of the probability shifts 
away from the rough surface, and reaches its maximal 
value for Nt — ► oo (practically, for Nt — 100). 

The eigenstates with small values of Ir emerge in the 
energy spectrum for any energy, thus resulting in a very 
slow convergence to the limit of ergodicity. For exam- 
ple, the number of eigenstates originated from m = 1 
is N m= i(E) = ^-y/E, therefore, the fraction of such 

eigenstates is given by %^ = f*-^ = y/^j^ 

where N(E) is the total number of states with energy 
less than E. One should stress that in the energy spec- 
trum these states appear regularily, due to the expression 
Eq.(5) with m = 1 and different n. 




FIG. 4. Two eigenstates that are strongly localized in 
a;— direction; (a) a — 390, (b) a = 407. The model parame- 
ters are the same as in Fig. 1. 

We would like to note that the type of localization in 
the channel space we discuss here, is different from that 
studied for circular billiards with a rough surface (see, for 
example, [14]). The point is that in our case classical dif- 
fusion in the transverse momentum space turns out to be 
very strong compared with quantum localazation effects 
[15]. In contrast to circular billiards, in our model with 
quasi-lD geometry the effects of a strong localization (in 
the channel space) are due to the existence of a contin- 
uous set of classical horisontal trajectories ("bouncing 
balls") which do not touch the rough boundary. As is 
shown in Ref. [16], these trajectories result in anomalous 
properties of conductance fluctuations for open waveg- 
uides of finite length. 

It should be stressed that for any high energy one can 
find eigenstates of a very different structure. To demon- 
strate this fact, in Fig. 4 we report two typical examples 
of strongly localized (in the coordinate x) eigenstates. 
These eigenstates are widely spanned in the unperturbed 



basis of H , with large values Ir ^> 1, and they corre- 
spond to large values of m <C 1. 

To conclude with, we have analysed the structure of 
eigenstates of a quasi-lD waveguide with a rough surface, 
paying main attentsion to their localization properties in 
the channel and coordinate rcprcscntaion. Different sets 
of strongly delocalized eigenstates (along the waveguide) 
have been found, that have a quite specific form in the 
transverse direction. We have develop the approcah that 
allows one to explain this form, using the transformation 
to new canonical variables. 

Another result is that eigenstates turn out to have very 
different localization properties and this difference can 
not be treated as a result of fluctuations only. Apart 
from the fluctuations, there are regular effects which are 
due to a strong influence of the geometry. Namely, it was 
found that the eigenstates originated from small values 
of to are strongly localized (delocalized) in the channel 
(coordinate) representation, and those associated with 
large m are strongly delocalized (localized) ones. This 
effect seems to be directly related to that found in Rcfs. 
[5,17] for open finite waveguides, where it was shown that 
characteristic scales for scattering are different for differ- 
ent channels. Thus, it seems questionable whether stan- 
dard statistical approaches based on completely random 
mathematical models, can adequately describe properties 
of eigenstates. 
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